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1.0  INTRODUCTION 

The  adaptive  wall  concept  is  a  systematic  procedure  for  eliminating  wall  interference. 
Ferri  and  Baronti  {Ref.  1)  and  Sears  (Ref.  2)  independently  originated  the  adaptive  wall 
concept  wherein  it  was  recognized  that  by  measuring  two  flow  variables  and  evaluating 
functional  relationships  for  unconfined  flow  conditions,  one  can  adjust  the  local  wall 
boundary  to  achieve  interference-free  flow.  The  procedure  is  iterative,  and  wall  boundary 
adjustments  are  made  until  unconfined  flow  is  achieved  at  a  measuring  surface  near  the  wall. 
The  analytical  investigations  by  Lo  and  Kraft  (Ref.  3)  and  Sears  (Ref.  4)  have  shown  the 
feasibility  of  the  concept.  Also,  a  numerical  simulation  by  Erickson  and  Nenni  (Ref.  5) 
established  the  concept,  and  recent  experimental  studies  (Refs.  6  through  9)  have  validated 
the  practicality  of  an  adaptive  wall  tunnel.  ' 

Since  the  adaptive  wall  procedure  is  iterative,  the  reduction  of  the  number  of  wall 
adjustments  to  achieve  unconfined  flow  would  save  on  tunnel  run  time  and  computer  time  in 
an  adaptive  wall  wind  tunnel.  This  has  been  the  impetus  for  the  study  of  accelerating 
convergence  and  the  development  of  one-step  convergence  formulae  to  evaluate  the  flow 
conditions.  The  one-step  formulae  significantly  reduce  the  number  of  iterations,  and  under 
certain  conditions  unconfined  flow  can  be  achieved  in  one  single  step. 

This  report  is  concerned  with  the  development  of  methods  for  use  in  two-dimensional 
and  axisymmetric  adaptive  wall  wind  tunnels.  Convergence  of  the  iterative  procedure  and 
derivation  of  the  one-step  convergence  formulae  are  presented  for  both  two-dimensional 
and  axisymmetric  flows  using  linearized  small  disturbance  theory  and  the  Fourier  transform 
technique. 

Numerical  examples  of  the  method  in  which  a  computer  program  serves  as  an  analog  of 
the  wind  tunnel  to  simulate  the  interior  flow  field  are  presented.  These  numerical 
demonstrations  have  shown  the  validity  of  the  adaptive  wall  procedure  and  the  power  of  the 
one-step  formulae  in  reducing  the  number  of  wall  iterations  in  subsonic  and  transonic  flows. 

2.0  ITERATIVE  PROCEDURE 

The  concept  of  the  adaptive  wall  wind  tunnel,  as  explained  in  Refs.  1  and  2,  states  that  to 
determine  whether  unconfined  flight  conditions  are  obtained  in  a  tunnel  with  any  model,  it 
is  necessary  and  sufficient  to  match  the  measured  flow  variables  at  a  selected  control 
surface,  say  S,  with  the  flow  variables  which  satisfy  the  unconfined  boundary  condition. 
Specifically,  it  is  necessary  to  measure  the  distributions  of  two  independent  flow  variables 
such  as  the  longitudinal  and  normal  perturbation  velocity  components  at  the  surface  S 
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(which  is  located  inside  the  tunnel  but  away  from  the  model).  One  of  the  measured  variables 
is  used  as  the  boundary  value  to  uniquely  specify  the  flow  field  exterior  to  S  in  the  presence 
of  the  condition  of  unconfined,  undisturbed  flow  of  a  uniform  stream  at  infinity.  Hence, 
since  the  two  measured  distributions  constitute  redundant  boundary  data  in  the  presence  of 
the  exterior  region  far-field  boundary  condition,  equivalence  at  S  of  the  measured  flow 
variables  interior  to  S  and  the  computed  flow  variables  exterior  to  S  constitutes  a  definition 
of  interference-free  flow  in  the  wind  tunnel.  Therefore,  by  comparing  the  exterior  region 
calculated  values  with  the  measured  values  of  the  same  quantities,  one  can  determine 
whether  or  not  unconfined  flow  conditions  exist. 

Unconfined  flow  conditions  can  be  achieved  if  provisions  are  available  for  adjusting  the 
wall  in  a  logical  manner.  A  basic  iterative  procedure  for  adjusting  the  walls  to  achieve 
unconfined  flow  is  presented  in  block  diagram  form  in  Fig.  1  and  applies  for  both  two-  and 
three-dimensional  flow  fields.  First  a  flow  field  is  established  in  the  tunnel,  and  the  velocity 
components  uT  and  vT  are  measured  at  the  given  control  surface  S.  The  exterior  unconfined 
region  is  then  evaluated  by  specifying  vE  =  vT  as  the  boundary  value  at  S.  If  the  distribution 
at  S  of  uE  determined  from  the  exterior  region  calculation  does  not  agree  with  uT,  then  the 
flow  is  still  affected  by  the  walls  and  the  walls  must  be  readjusted.  The  iteration  continues 
until  uE  and  uT  agree.  Then  the  flow  about  the  model  in  the  tunnel  is  unconfined. 

To  accelerate  convergence,  a  relaxation  factor,  k,  is  introduced  in  the  iterative  procedure 
as 

uT(n+1)  =  kuE(n)+  (1  -k)uT(n)  (1) 

at  S.  The  tunnel  wall  is  adjusted  according  to  Eq.  (1).  The  proper  selection  for  the  value  of  k 
to  speed  up  iteration  convergence  will  be  discussed  later. 

Alternately,  uE  =  uy  can  be  specified  as  the  boundary  condition  on  S  in  the  exterior 
region,  and  vE  can  be  compared  with  Vy  to  determine  whether  unconfined  flow  exists  in  the 
tunnel.  More  generally,  any  two  independent  measurable  flow  variables  can  be  used  in  the 
adaptive  wall  procedure. 

3.0  CONVERGENCE  OF  THE  ITERATIVE  PROCEDURE 

To  provide  an  analytical  proof  of  convergence  of  the  adaptive  wall  iterative  procedure, 
an.  interior  tunnel  simulation  is  required  in  addition  to  an  exterior  flow  functional 
relationship. 
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3.1  TWO-DIMENSIONAL  TUNNEL 


3.1.1  Exterior  Flow-Field  Relationships 


As  indicated  previously,  an  exterior  flow-field  relationship  is  required  to  demonstrate  the 
adaptive  wall  concept.  The  flow  field  is  described  by  the  linearized  Prandtl-Glauert 
equation,  and  the  boundary  value  problem  in  the  exterior,  unconfined  flow  region  is  shown 
in  Fig.  2.  The  boundary-value  problem  can  be  solved  by  the  Fourier  transform  technique 
using  vE(x,  ±h)  as  the  boundary  value  at  y  =  ±h  along  with  the  condition  of  unconfined 
flow  at  infinity.  The  corresponding  unconfined  flow  distribution  Ue(x,  ±  h)  at  the  measuring 
plane  S  is  related  to  the  vE(x,  ±  h)  in  the  transformed  plane  as 


—  1  p  — 

uE(p,±M  =  -  —  vE(p,±h) 

P  IpI 


where  the  barred  quantities  are  the  Fourier-transformed  variables  defined  by 

g(p,±h)  =■  (277rl/2  /  g(x,±h)  elpK  dx 


(2) 


/ 


with  p  being  the  Fourier  transform  parameter. 

Alternately,  if  uE  is  the  boundary  condition  at  S  for  the  exterior  region,  then 

VE  =  ~'lP  TT  ue(R*m 

■  IpI 


(3) 


Equations  (2)  and  (3)  can  be  inverted  into  the  physical  plane  as 


which  are  the  familiar  exterior  flow  functional  relationships  presented  in  Ref.  2. 
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3.1.2  Interior  Flow-Field  Simulation  .  .  / 

Simulation  of  the  interior  flow  field  is  required  strictly  for  the  theoretical  proof  of 
convergence.  It  should  be  emphasized  that  the  power  of  the  adaptive  wall  concept  is  that  it 
never  requires  the  calculation  of  the  interior  flow  field  when  applied  to  a  real  wind  tunnel. 
The  measured  tunnel  flow  is  simulated  here  by  the  analytical  model  of  the  flow  in  the  tunnel 
described  by  the  boundary- value  problem  shown  in  Fig.  2. 

Based  on  thin  airfoil  theory,  the  thickness  and  camber  (including  incidence)  of  an  airfoil 
may  be  treated  as  symmetric  and  asymmetric  problems,  respectively,  if  the  control  surfaces 
are  selected  symmetrically  about  the  airfoil.  The  boundary-value  problem  outlined  in  Fig.  2 
can  be  divided  into  symmetrical  and  asymmetrical  problems  as  given  in  Fig.  3.  Thus,  the 
velocity  components  at  the  control  surface  can  also  be  divided  into  symmetrical  and 
asymmetrical  components  as 

uT(±h)  =  Up  (±h)  +  uT  t±li)  (6) 

s  A 

and 

vp  (ill)  =  vT  (+h)  +  VT  (±h).  (7) 

8  A 

Then  the  symmetric  and  asymmetric  perturbation  velocities  can  be  written  in  terms  of  the 
velocity  components  at  the  control  surfaces  as 


Up  (±h)  = 

‘  s 

1  j 

2  1 

[vip  (h)  +  Up  ( — h)  J 

(8a) 

vT  (+h)  = 

1  5 

'1 

+  — 

■  2 

j\'p{h)  -  vT(— h)] 

(8b) 

uT  (±h)  = 
n 

1 

t  2 

[uT  (h)  -  up(-h)] 

(9a) 

vt  (+b)  = 

A 

1 

2 

[  vp  (li)  +  v p  (— h ) J 

(9b) 

For  the  symmetric  problem,  the  transformed  normal  velocity,  vTs,  in  the  tunnel  can  be 
expressed  as  a  function  of  the  model  geometry  and  an  arbitrary  distribution  of  the 
longitudinal  velocity,  ups ,  as 

vp  (p,h)  =  i/3  —  tanh  (|pj/3h)  up  (p^h)  -  sech  (p/3h)  F(p)  (10) 

S  |  P I  S 

/ 
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where 


F(p) 


(a,)'1'2 


(x)  eipx.  dx  ; 
■  \ 


The  model  thickness  distribution,  F(x),  can  be  interpreted  as  a  “potential  equivalent” 
airfoil  profile  including  viscous  effects.  It  is  reasonable  to  assume  that  the  same  equivalent 
profile  exists  in  an  unconfined  flow  condition.  The  solution  for  unconfined  flow  is  readily 
obtained  in  the  transformed  plane  (Ref.  3)  as 

(P.h)  =  '--^r  e'|p|^p>  (,1) 

*  p  IpI 

vM  «p,h)  =  e'lpl  £h  F(P)  (12) 


For  the  asymmetric  problem,  the  normal  velocity  in  the  tunnel  can  be  expressed  as  a 
function  of  the  model  camber  and  incidence  and  as  an  arbitrary  distribution  of  the 
longitudinal  velocity;  thus, 

VT  (p,h)  =  i|8  coth  (|p||8h)  UT  (p,h)  -  csch  (| p|/3h)  y  (p)  ,  (*3) 


I P  1 


2  IpI 


where 


OG- 

'  (p)  =  (2tt)'^  £  y(x)elpxdx 


The  distribution  of  vortices,  7(x),  is  used  to  represent  the  model  camber  and  incidence.  With 
the  assumption  of  vortex  strength  remaining  the  same  for  confined  and  unconfined  flow 
conditions,  the  solution  for  unconfined  flow  is  expressed  as 


uoc  fpd1)  =  —  e  ^pl^y(p) 


(14) 


-  .(  1\  P  —  ( p  I  jQ h 

V"A(P'h)  -  .Th’ 


y(p) 


05) 


3.1.3  Proof  of  Convergence  ' 

The  convergence  of  the  iteration  procedure  will  be  proved  for  the  symmetric  problem 
first  in  the  transformed  plane  since  the  functional  relationships  and  tunnel  simulation  are  in 
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simple  algebraic  form  in  that  plane.  The  iteration  is  initiated  with  an  unspecified  tunnel  wall 
r  ,  ,  (o)  (o) 

configuration  which  provides  the  arbitrary  initial  value  of  flow  quantities,  Ujs  and  vTs,  or  in 
the  transformed  plane,  uE  and  v\°\  at  the  measuring  plane  y  =  h.  Following  the  iterative 
scheme  shown  in  Fig.  1,  one  can  obtain  the  nth  iterative  values  for  ues  and  Uts  at 
y  =  h  by  appropriately  utilizing  Eqs.  (2),  (3),  and  (10)  to  yield 


and 


(n) 

UE  (p,h) 


p  _>>  i 

—  tanh  (p/9h)u™  (p,h)  +  -A,n  =  0 

IpI  8  0 


—  tanh  (pj8h)  G*n^  +-  —  A, 

IpI  p 


1,  2, 


u  tp,h)  =  w  u  E  (p,h)  +  (1  -  a>)  u  j  (p,h) 
r»  s  8 


=  G(n>,n  =  1,  2, 


.(17) 


where  lu  is  a  relaxation  factor  defined  in  the  transformed  plane  to  accelerate  convergence  and 

(i-n"-1) 


=  Qii-i  GU)+  lw 


(3  (l—O) 


A ,  n  =  2,3, 


where 


(18) 


.  (o)  Jo) 

G*-1^  =  —  -vT  (p,h)  +  (1  —  <u)  u-p  (p,h) 

P  'pi  S  s 

fl  =l—iu  +  tanh  (|p  |  $h)J 


(19) 

(20) 


A  =  JL-  sech  (p/3h)  F  (p). 

IpI 


(21) 


It  should  be  noted  that  fl2  <  1  for  0  <  oj  <  1.  As  the  iteration  proceeds,  the  effect  of  the 
initial  tunnel  conditions  (represented  by  GO))  diminishes.  In  the  limit  as  n  —  «  it  is  readily 
established  that 

fim  c(n)  =  ..  i<oA  =  L  e-|pl^hF( p)  (22) 

n  *  “  |8(i-Q)  /8  |p| 
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which  is  equal  to  u^p.h)  as  indicated  in  Eq.  (11).  Consequently,  in  the, limit  n  —  ». 


(n) 

UE  (p,K) 

S 


£im 

n  » 


Jn.) 

Uj  (p,h)  = 

a 


eim  G(n)  =  u  (P,h)' 

l  ->  no  c*o  r 7 

s 


(23) 


and,  hence,  convergence  of  the  adaptive  wall  method  to  unconfined  flow  is  established 
independent  of  the  initial  tunnel  conditions  for  nonlifting  airfoils  in  subsonic  flow.  It  is 
important  to  note  that  the  iteration  procedure  converges  independently  of  the  relaxation 
factor  for  values  0  <  u  <  1 .  It  should  also  be  noted  that  convergence  can  be  demonstrated 
in  a  similar  manner  for  the  alternate  procedure  of  specifiying  v-ps  instead  of  Uys  at  the  control 
surface. 


Similarly,  the  asymmetric  problem  can  be  shown  by  using  Eqs.  (2),  (3),  and  (13) 
following  the  scheme  shown  in  Fig.  1.  Since  the  problem  is  linear,  the  general  problem 
including  the  symmetric  and  asymmetric  portion  is  also  convergent  independent  of  the  initial 
tunnel  conditions. 

3.2  AXISYMMETRIC  TUNNEL 


3.2.1  Exterior  Flow-Field  Relationships 


The  flow  field  in  the  region  exterior  to  and  including  a  circular  measuring  surface  at 
r  =  R  is  described  by  the  linear  potential  flow  equation  in  axisymmetric  coordinates  and  is 
shown  in  Fig.  4.  With  the  boundary  condition  specified  at  r  =  R  as  vE(x,R)  and  at  infinity  as 
unconfined  flow,  the  potential  flow  equation  can  be  solved  by  Fourier  transforms  for  the 
exterior  region.  The  corresponding  flow  distribution,  uE(x,R),  at  the  measuring  surface 
(r  =  R)  is  related  to  vE(x,R)  in  the  transformed  plane  by 


UE 


p  MpI0r) 
71  Kt  (!  P  iySR) 


(24) 


where  K0  and  Ki  are  modified  Bessel  functions  of  zero  and  first  order. 
If  uE  is  specified  as  the  boundary  condition  at  r  =  .R,  then 


vE(p,R)  =  - 


p  Kx  (|p|jBB> 

H  MpP 


uE  (p,R) 


(25) 
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The  unconfined  flow  field  is  determined  about  a  body  of  revolution  with  its  axis  at  r  =  0. 
Using  the  Fourier  transform  technique,  one  can  derive  the  streamwise  and  normal  velocity 
components  at  the  measuring  surface  as 

uM  (P,R)  =  ip  Ko(|p  |jQR)  Q (p)  ^26) 

and 


(p,R>  -  |p |  K1(|p|0R)Q(p)  (27) 


where 


Q(P)  =  (2n)-l/> 


U)  e1?*  dx 


3.2.2  Interior  Flow-Field  Simulation 

In  order  to  verify  convergence  of  the  iterative  procedure,  it  is  necessary  to 
mathematically  simulate  the  interior  flow  field  as  shown  in  Fig.  4.  If  the  linear  potential  flow 
equation  in  the  transformed  plane  is  used  with  boundary  conditions  specified  as  vx(p,R)  and 
the  model  geometry  as  used  in  Eq.  (26),  the  solution  to  the  boundary-value  problem  for  the 
interior  region  can  be  derived  as 


u  r  (p,  R) 


_j_  _P_  _P_  Q  (p) 

P  |p|  MIpI0R)  Vt  P’  +  IpI  r  MpI^ 


(28) 


where  I0  and  I|  are  modified  Bessel  functions.  If  uT(p,R)  is  specified,  the  solution  becomes 


vT 


t^IpIiSr)  _  Q(p) 

-  Hrr(p,R}+  - 

l0(|p|/3R>  RIo(|p|/3R) 


(29) 


3.2.3  Proof  of  Convergence 

In  the  transformed  plane,  convergence  of  the  iterative  procedure  shown  in  Fig.  1  can  be 
demonstrated.  Using  Eqs.  (24),  (25),  (28),  and  (29),  the  nth  iterative  values  of  uE(p,R)  and 
ur(p,R)  are  determined  from  given  initial  values,  ux}  (p,R)  and  v^  (p,R),  to  be 
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K(|p|/3R)  l,(|pli8R)  _>).  i 

- ;  UT  (p,R)  +  -  A»  n 

k^IpI/Sr)  ib(|p|^R)  T  0 


=  0 


(n) 

uE(p,R)  = 


Ko  CIpIjSR)  ^CIp^R)  ^  ^  i  a 

KI  C|p|(8R)  io(|p|j8R)  +  P 


(30) 


n  =  1,  2, 


and 


(n)  _ (n- 1 ) 

uT  (p,R)  =  a>  uE  (p,R)  +  (1  —  o>)  uT  (p»R) 


=  G(n\  n  =  1,  2, 


(31) 


where  u  is  a  relaxation  factor  defined  in  the  transformed  plane  and  is  selected  to  accelerate 
convergence.  The  recurrence  relationship  of  the  iterative  scheme  is 


)n-  1 


where 


G(n)  =  Qn-1Gm  +  i  —  1 . .  A,  n  =  2,3, 

j8  l  -  ft 


(o)  _(o) 

;(1)  =  i  —  vj  (p,R)  +  (1  —  tu)  uT  (p,R) 


(32) 


(33) 


Q 


1  - 


1  + 


Ko  (|p|jSR)  ijdpl^R)! 
K^IpI/SR)  l„(|p|j8R)_ 


(34) 


!  Ko(|p|j8R)  Q  (p) 
r  I0(|pl/3R)  k,  (|p l/SR) 


(35) 


It  can  be  shown  that  O2  <  1  for  0  <  w  <  1.  With  this  restriction  on  qj,  the  effects  of  the 
initial  conditions  diminish  as  the  iteration  proceeds.  In  the  limit  as  n  —  »,  the  recurrence 
relationship  becomes 

ipKo(|p|0R)Q(p) 


Pim  G(n) 


>\  1 


j8  1  -  Q 
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which  is  equal  to  uB  (p,R)  as  indicated  in  Eq.  (26).  Also,  in  the  limit  as  n  — 

Jn)  Jn) 

lim  uE  (p,  R)  =  lim  u  -p  (p,R)  =  lim  G^‘  =  u  (p,  R) 

and,  hence,  convergence  for  the  axisymmetric  iterative  scheme  has  been  established  which  is 
independent  of  the  initial  conditions. 

4.0  ONE-STEP  CONVERGENCE  FORMULAE 

A  technique  to  speed  up  the  convergence  of  the  iteration  process  is  important,  as 
indicated  before.  For  subsonic  flow,  a  formula  can  be  derived  for  determining  conditions  of 
unconfined  flow  in  one  single  step.  Consequently,  for  an  adaptive  wall  wind  tunnel  it  is 
conceivable  that  unconfined  flow  could  be  attained  in  a  single  adjustment  of  the  tunnel 
boundary. 

4.1  TWO-DIMENSIONAL  TUNNEL 

The  conditions  for  unconfined  flow  in  a  tunnel  can  be  achieved  by  properly  choosing  an 
optimal  relaxation  factor.  From  Eq.  (20),  the  optimal  relaxation  factor,  u>,  can  be  obtained 
by  setting  ft  =  0  as 

l 

OJ)t  1  +  tanh  (] p  |(3h)  (36) 


which  is  in  the  transformed  plane.  One-step  convergence  formulae  are  derived  by 
substituting  wopt  of  Eq.  (36)  into  Eq.  (17)  as  given  in  Ref.  3  for  the  symmetric  case. 


Alternately,  one  can  obtain  one-step  formulae  to  relate  the  unconfined  flow  variables 
with  tunnel  variables  through  the  model  profiles.  Specifically,  the  elimination  of  F(p) 
between  Eqs.  (10),  (11),  and  (12)  yields  the  one-step  formulae  for  the  symmetric  problem  as 


(p,h) 


,-|pl£h 


_  1  p 

sinh  (|p  | jQh)  u^,  (p,h)n - - —  cosh  (p/3h) 

■  0  Pi 


v-p  (p,h) 

•S 


(37) 


and 


v  (p,h)  = 


■l£h 


sinh  ( [ p  |j8h)  uT  (p,h)  —  cosh(p/3h)  vT  (p,h) 


(38) 


s 
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Similarly,  one  combines  Eqs.  (13),  (14),  and  (15)  to  eliminate  the  distribution  of  vortices, 
7(p),  and  obtains  the  one-step  formulae  directly  related  to  the  tunnel  flow  variables,  uTa  and 
vTa  ,  at  the  upper  control  surfaces  as 


u  (p,h)  =  e“lp1^h  [cosh  (p/9h)  ux  (p.h)  +  ~  ~  sinh  (|p  |jBh)  vT  (p,h) 

“a  |_  a  p  Ip  I  A 


(39) 


and 


v,*,  (p,h)  =  -e‘ 


=  _e-|p|$h  ifS  — — -  cosh  (p$h)  uT  (p,  h)  -  sinh  (jp  |/3h)  vT  (p,  h)|  (40) 

Ip! 


,h) 


For  a  thin  airfoil,  the  desired  one-step  formulae  are  obtained  by  superposition  of 
thickness  and  camber  effects  from  Eqs.  (37),  (38),  (39),  and  (40)  with  the  aid  of  Eqs.  (8)  and 
(9). 

U  <*,  (h>  =  u»o  W 

a.  A 


.(h) 


.  1  e-|p|(3h 


|p|£h  uT(h)-e~|pl^h  uT(-h) 


+  i-L  e|pi0h  v  (h)  _  i  JL  e-|pl^7T(-h)1 
£  I p |  T  J8  IpI  T  j 

=  _  i..  e  |p|(Sh  r_  e~!  P  l^h  vT(h)  +  Jp'/31*  7T(-h) 

+  i/3  —  elP l^h  uT (K)  -  i/3  e"|p|^h  uT(-.h)l 

1 p I  ■  IpI  J 


(41) 


(42) 


Similarly,  u^-h),  vM(h)  at  the  lower  control  surface  can  also  be  obtained.  The  above 
pair  of  equations  can  be  inverted  into  the  physical  plane  as  the  following  expressions: 

f  X  /3h  r  *T^'+h) 

u  U,  ±h)  =  —  uT  (x,  ±h)  — —  J  —TZ - — 


k(£-X> 


2iT/S 


(f,  ±h)  , 

-  a£  — 

f-X  277$ 


f 


v T  (£,  T  h) 

k(£-*) 


(f-x)  ^ 


(43) 
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and 


where 


Voe  (x,ih)  =  -  vT  (x,±h)  -  —  • 


/ 


°°  VT  (£,  T  h) 


k(£-x). 


it 


+ _ 

2  n 


8  ?  uT(£±h)  8  i°  uT^’Th) 

/  f  -77T-T-. 


k(f-x) 


k(f-  x)  =  (2/3h) 2  +  (f-x)2 


(44) 


and  uw,  u-f,  v^,  and  vT  are  the  perturbation  flow  variables.  The  preceding  one-step 
convergence  formulae  give  the  unconfined  flow  variables,  uOT(x,  ±h)  and  vOT(x,  ±h)  at  the 
control  surface,  when  the  flow  variables,  uT(x,  +  h)  and  vT(x,±h),  are  in  hand  for  any 
arbitrary  model. 

4.2  AXISYMMETRIC  TUNNEL 

As  in  the  two-dimensional  case,  the  one-step  formulae  for  accelerating  convergence  can 
be  derived  by  relating  tunnel  variables  in  Eqs.  (28)  and  (29)  with  the  unconfined  variables  in 
Eqs.  (26)  and  (27)  and  eliminating  the  model  geometry  Q(p)  to  obtain 


u„c(p.R)  =  ip(3R  ko  (|pl£R) 


«0  (|P|i8R) 


vT(p,  R)  -  i  I1(|p|/3R)uT(p,R)J 


(45) 


and 


v«,  (p.  R)  =  IpIjSRMpI/SR)  I0(|p|/3R)  vT(p,R)  -  ijS  — Ijdp'l/SR)  uT(p,R) 


(46) 


An  approximate  expression  in  the  physical  plane  is  derived  for  Eq.  (45)  in  Appendix  A.  The 
resulting  expression  is 


(x,  R)  =  -  j^uT(x,  R)  ■+  uE(x,R)j 


M 

£  / 

1=1  ^D,° 


vT  (£  R)  -{£-*) 

- -  df 

(£-*)2  +  (mj0H/2)2 


1 

rr 


M 

L 

m=  1 


c 

m 


oa  (£,  R)  •  (mjSR/2) 

/  - -  it' 

-Z  (^-x)2  +  (mj8R/2)2 


(47) 
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\ 


t 

where  Ug  is  the  exterior  flow  functional  relationship  equation,  Eq.  (24), 


u£(x,  R)  =  - 


1 

rr/3 


/ 


— OC 


vT  (£  R) 

f-x 


M 


I 


vT  {f,  R)  ■  (<f  —  x) 

— — -  d£ 

(f-x)2+  (m/SR/2)2 


where  am,  bm,  and  cm  are  coefficients  of  exponential  curve  fits  to  Bessel  functions  defined  in 
Appendix  A, 

An  expression  for  the  other  unconfined  variables  vw(x,R)  has  not  been  inverted  into  the 
physical  plane.  In  the  following  numerical  examples,  this  variable  is  not  required  directly  but 
can  be  determined  from  the  finite-difference  computation  on  the  next  iteration  using 
u^fx.R)  as  the  boundary  condition  for  the  interior  flow  field. 

5.0  NUMERICAL  EXAMPLES 
5.1  TWO-DIMENSIONAL  TUNNEL 


5.1.1  Subcrilical  Flow 

An  NACA  0012  airfoil  at  zero  angle  of  attack  and  M  =  0.6  is  chosen  as  the  first  example 
to  illustrate  the  convergence  of  the  adaptive  wall  iterative  procedure.  For  simplicity,  the 
open-jet  boundary  is  selected  at  the  control  surface  symmetrically;  i.e.,  u?(x,±h)  =  0.  The 
interior  flow  field  is  computed  by  the  transonic  airfoil  program  TSFOIL  (Ref.  10)  to  obtain 
v($x,  ±h)  as  shown  in  Fig.  5  by  a  dashed  line.  Notice  that  the  normal  perturbation  velocity 
component  is  converted  to  degrees  to  make  it  consistent  with  units  of  flow  angle. 

With  the  distributions  of  Uj(x,  ±  h),  v?\x,  ±  h),  the  basic  iterative  procedure  outlined  in 
Fig.  1  is  carried  out  using  Eq.  (1)  with  k  =  1.0.  The  boundary  condition  is  updated  on  each 
successive  iteration  by  the  exterior  relationship  in  Eq.  (4)  and  convergence  to  the  unconfined 
distributions  is  achieved  by  the  fifth  iteration.  The  results  for  the  velocity  components  at  the 
control  surface  are  presented  in  Fig.  5  and  for  the  airfoil  surface  pressure  distribution  in 
Fig.  6.  1 

Using  one-step  formulae,  Eqs.  (43)  and  (44),  with  the  same  initial  distributions  of 
Uy!(x,±h)  and  Vj’(x, ±h),  one  can  determine  the  u00(x, ±h)  and  v^fx^h)  at  the  control 
surface  in  a  single  step.  These  distributions  are  in  good  agreement  with  a  separate 
computation  with  TSFOIL  for  unconfined  flow  as  shown  in  Fig.  7,  which  demonstrates  the 
capability  of  the  one-step  formulae.  Furthermore,  with  the  boundary  conditions  specified  as 
Uoo(±h)  or  v^fih)  in  the  tunnel  region,  the  airfoil  surface  pressure  distribution  is  obtained 
as  showm  in  Fig.  8,  which  shows  excellent  agreement  with  the  unconfined  solution. 
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The  second  example  is  a  lifting  case  with  a  NACA  001 2  airfoil  at  a  =  I  deg  and  M  =  0.6. 
Starting  from  a  closed  tunnel,  i.e.,  Vj[x,  ±h)  -  0  at  the  control  surface,  the  distribution  of 
Uj(x,  ±h)  using  TSFOIL  is  presented  in  Fig.  9.  With  the  distribution  of  Uy|x,  ±h), 
v^|x,±h),  the  one-step  equation  determines  u^fx^h),  v^fx, ±h)  which  agrees  with  a 
separate  TSFOIL  computation  as  shown  in  Fig.  9.  Similarly,  the  airfoil  pressure  distribution 
illustrates  the  excellent  agreement  with  unconfined  flow  results  as  shown  in  Fig.  10.  ■ .  . 

5.1.2  Supercritical  Flow 

For  subsonic  Mach  numbers  with  locally  supercritical  flow,  the  one-step  formulae  can  be 
used  in  the  iteration  procedure  to  accelerate  convergence.  An  NACA  0012  airfoil  at 
M  =  0.8,  a  -  0  yields  results  as  shown  in  Fig.  11  with  only  one  additional  iteration  to 
achieve  convergence.  For  lifting  conditions  M  =  0.72  and  a  —  1  deg,  convergence  is  also 
achieved  by  the  second  iteration  as  shown  in  Fig.  12. 

5.2  AXISYMMETRIC  TUNNEL  -  "  ■ 

Two  numerical  examples  chosen  to  demonstrate  convergence  of  the  iterative  procedure 
for  axisymmetric  flow'  have  utilized  the  one-step  expression  of  Eq.  (47).  In  both  examples  a 
transonic  small  disturbance  program  is  used  to  simulate  the  interior  tunnel  flow  field  over  a 
parabolic  arc  body  of  revolution  with  fineness  ratio  of  ten.  The  boundary  condition,  uy(R) 
at  the  control  surface,  is  varied  on  each  successive  iteration  based  on  the  calculations  from 
Eq.  (47).  In  the  first  example,  the  flow  is  subcritical  throughout,  and  in  the  second  example 
the  flow  field  is  supercritical  in  the  vicinity  of  the  model  but  subcritical  at  the  control 
surface. 

(o) 

In  the  first  example,  Mach  number  is  0.9,  R/f  is  0.5,  and  the  initial  condition,  uT(x,R), 
represents  an  open  jet  boundary  condition.  The  flow'  field  is  subcritical  and  remains 
subcritical  when  the  boundary  conditions  are  updated  from  Eq.  (47).  Therefore, 
convergence  to  the  unconfined  flow  conditions  using  the  one-step  formula  should  be 
achieved  on  the  first  iteration.  The  results  in  Fig.  13  verify  that  convergence  has  been 
reached  on  the  first  iteration  for  the  control  surface  velocity  component  distributions  and 
the  body  surface  pressure  distribution. 

In  the  second  example,  Mach  number  is  0.975,  R/fis  1.0,  and  the  initial  condition, 
ut(x,R),  is  selected  as  an  open  jet  boundary  condition.  Because  of  the  increased  Mach 
number,  the  flow  field  becomes  supercritical  but  remains  subcritical  at  the  control  surface. 
Under  these  conditions,  the  one-step  convergence  formulae  accelerate  convergence. 
Convergence  is  achieved  on  the  third  iteration  as  indicated  in  Fig.  14. 
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6.0  CONCLUDING  REMARKS 

The  analytic  proof  of  convergence  of  the  adaptive  wall  method  using  the  subsonic  small 
perturbation  theory  has  established  the  validity  of  the  concept  of  the  adaptive  wall  wind 
tunnel.  The  development  of  one-step  convergence  formulae  to  increase  the  convergence 
speed  of  the  iterative  procedure  has  provided  the  practicality  of  the  adaptive  wall  method. 
Numerical  simulations  have  demonstrated  the  effectiveness  of  one-step  formulae  in 
accelerating  convergence  for  both  subsonic  and  transonic  flows. 
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Figure  1.  Block  diagram  of  the  iterative  procedure  of  the  adaptive  wall  concept. 
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Figure  2.  Two-dimensional  boundary-value  problem  of  exterior 
unconfined  region  and  interior  tunnel  region. 
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a.  Thickness  problem  b.  Lifting  problem 

Figure  3.  Mathematical  model  for  thickness  and  lifting  problems. 
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Figure  4.  Axisymmetric  boundary-value  problem  of  exterior 
unconfined  region  and  interior  tunnel  region. 


AEDC-TR-79-55 


AEDC-TR-79-55 


o  UNCONFINED  SOLUTION 

-  CLOSED  TUNNEL 

-  ONE- STEP  FORMULAE 


x/c 


Figure  9.  Velocity  component  distributions  from  one-step  formulae 
at  the  control  surface  for  NACA  0012  airfoil,  h/c  =  1.0, 
M  =  0.6,  a  =  1.0; 
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Figure  10.  Airfoil  surface  pressure  distributions  based  on  one-step 
results  for  NACA  0012  airfoil,  h/c  =  1.0,  M  =  0.6, 
a  =  1.0. 
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Figure  13.  Distributions  of  velocity  components  at  the  control  surface  (R/fi  =  0.5) 

and  body  surface  pressure  on  a  parabolic  arc  body  of  revolution.  C/d  =  10, 
M  =  0.9,  a  -  0. 
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APPENDIX  A 


INVERSION  OF  THE  AXISYMMETRIC  ONE-STEP  CONVERGENCE  FORMULAE 

Unlike  the  case  of  two-dimensional  flow,  an  analytic  inversion  of  the  one-step  con¬ 
vergence  formulae  for  axisymmetric  flow,  Eq.  (45),  into  the  physical  plane  is  not  available. 
An  approximate  inversion  is  derived  as  shown  below.  Using  the  exterior  flow  functional 
relationship,  Eq.  (24), 

.  Ko(|p|/8R) 

■  '  "e1"-11*  -  ^  Vt1p’h\ 

the  one-step  equation  in  the  transformal  plane,  Eq.  (45),  may  be  rewritten  in  the  form 

,  r_  _  qip  Ti  1- 

(P-  “  2  [uT(P’R)  +  uE(p’R)J+  ~  [j  “  (A)  (A)J  VT  (p’  R) 

_ri_XK0(\)l,  (A)J  uT(p,  R) 

(All  Bessel  functions  1,,  Ko,  K,  have  argument  X  =  |  p  |j0R.) 

This  can  be  inverted  into  the  physical  plane  as 
uM(x,R)  =  j|uT  (x,R)  +  uE  (x,R)j 

;  j 


(£  -  x) 


(£_  x)*  +  <mfSR/2 V 


R)  df 


where 


ur(x,R) 


4E-  /' 


m/3R/2 

(f  -x)2  +  (m/3R/2)2 


/oo  A  o 

E-  / 


ut(£R>  d£ 


(£_*) 


(£_x)2  +  (m/SR/2)2 


vx(£R)  df 


The  am,  bm,  and  cm  are,  respectively,  the  coefficients  of  the  exponential  series  approximation 
of  the  following  Bessel  functions: 

M 

K  (X)  —  K  {A) 

D  1 


K1(X) 


E 

m  =  1 


— mA/2 
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KoU) 
K ,  (A) 
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2  ° 


>w]  ■  s 

m=  1 


b  e 


— mA/2 


t-AK0U)  Ij(A)  =  >  c,  »-mA/2 


s 


m=l 


The  coefficients  am,  bm,  cm,  are  determined  by  truncating  the  curve  fits  at  M  -  5.  The 
value  of  these  coefficients  used  in  the  numerical  demonstrations  in  Section  5  are  as  follows: 

am  =  1-1.748,  9.557, -26.203,  31.026, -13.6241  -if; 
bm  =  jO.661, -1.979,  13.078, -0.609, -1.137} 

cm  =  10.835,  -3.379,  8.115,  -7.387,  2.323} 
m  =  1,  2,  .  .  ,  5 
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NOMENCLATURE 

Coefficients  of  exponential  curve  fits  (Appendix  A) 

Small  disturbance  pressure  coefficient,  -2u 
Airfoil  chord 

Maximum  diameter  of  axisymmetric  body 
2-D  model  profile 

Defined  in  Eq.  (17)  for  2-D  flow  and  Eq,  (31)  for  axisymmetric  flow 
Location  of  2-D  control  surface  (Fig.  2) 

Relaxation  factor  in  physical  plane 

Axisymmetric  body  length 

Fourier  transform  parameter 

Axisymmetric  model  profile 

Location  of  axisymmetric  control  surface  (Fig.  4) 

Control  surface 

Perturbation  velocity  components  parallel  and  normal  to  free  stream, 
respectively  (normalized  by  free-stream  velocity) 

Axisymmetric  coordinates  parallel  and  normal  to  free  stream,  respectively 
2-D  coordinates  parallel  and  normal  to  free  stream,  respectively 
Airfoil  incidence,  deg 
Vl  -  M2 

2-D  model  camber  and  incidence 

Defined  in  Eq.  (21)  for  2-D  flow  and  Eq.  (35)  for  axisymmetric  flow 
Perturbation  velocity  potential  normalized  by  free-stream  velocity 
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Q  Defined  in  Eq.  (20)  for  2-D  flow  and  Eq.  (34)  for  axisymmetric  flow 

u>  Relaxation  factor  in  transform  plane 

Subscripts 

A  Asymmetric  component 

E  Exterior  region 

m  Summation  index 

s  Symmetric  component 

T  Tunnel  interior  region 

<*>  Free-air  condition 

Superscripts  . 

(n)  Iteration  number 

*  Critical  value 

Notations 

|  |  Absolute  value 

(  )  Quantity  in  transform  plane 

£  Summation 
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